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Conservation of canonical angular momentum shows that charged particles are typically con- 
strained to stay within a poloidal Larmor radius of a poloidal magnetic flux surface. However, 
more detailed consideration shows that particles with a critical charge to mass ratio can have 
zero canonical angular momentum and so be both immune from centrifugal force and not con- 
strained to stay in the vicinity of a specific flux surface. Suitably charged dust grains can have 
zero canonical angular momentum and in the presence of a gravitational field will spiral inwards 
across poloidal magnetic surfaces toward the central object and accumulate. This accumulation 

results in a gravitationally-driven dynamo, i.e., a mechanism for converting gravitational potential 
i— i 

■ energy into a battery-like electric power source. 
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I. INTRODUCTION 



To the best of the author's knowledge, no plasma has ever been observed to make a Kepler 
orbit around the Sun, a planet, or a moon. For example, the Earth's magnetotail does not 
make Kepler orbits around the Earth nor does the solar wind make a Kepler orbit around the 
Sun. This is puzzling because in the astrophysical literature plasmas are routinely presumed 
to make near-Kepler orbits in the presence of the gravitational field of a central object, e.g., 
accretion disks are modeled as MHD plasmas in a near-Kepler orbit around a central object. 
The qualifier 'near' is used because it is conventionally presumed for gases and plasmas that 
a radial gradient of an isotropic pressure exists and that this pressure gradient provides 
a modest outward force which slightly reduces the amount of centrifugal force required to 
balance the inward force of gravity and so achieve a stable circular orbit (e.g., see Ref. [l| 
for a discussion of the implications of this effect in the context of a system consisting of gas 
and solid particles). 

If Kepler orbiting plasmas are so ubiquitous in astrophysics, then why is there not even 
a single example in the great variety of plasma/gravitational situations in our own solar 
system? A Kepler orbit is essentially a property of a single point particle - it is not a 
property of a collection of independent point particles. For example, Earth and Mars are 
individually in Kepler orbits around the Sun, but the center of mass of Earth and Mars is 
not in a Kepler orbit around the Sun. If one replaced Earth and Mars by some statistically 
large number of point particles, then each could be in its own Kepler orbit around the Sun, 
but the center of mass of this configuration would not be in a Kepler orbit. 

A collisional gas in a container could be in a Kepler orbit around the Sun because the 
walls of the container bind the particles to stay within a fixed distance of the center of 
mass so that the entire system can be considered as a point particle located at the center 
of mass. The transition from single point particle behavior to the behavior of a group of 
particles can be seen by considering a moon in close orbit around a massive planet. The 
tidal forces resulting from the gradient of the planet's gravitational force can be so strong 
as to overcome the binding forces and fracture the moon. The fragments would then follow 
distinct individual Kepler orbits and the center of mass of these orbits would not follow 
a Kepler orbit. Similarly, in order for a gas (or plasma) to have its center of mass follow 
a Kepler orbit there would have to be some binding force, such as would be provided by 
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container walls, that would prevent the particle constituents of the gas (or plasma) from 
separating onto distinct Kepler orbits. Collisions might provide binding at the interior of 
a gas cloud, but not at the periphery because particles at the periphery moving away from 
the center would not encounter other particles with which they could collide. The periphery 
would simply expand into vacuum if there is no wall to prevent this expansion and the 
particles constituting the periphery would make Kepler orbits substantially different from 
the Kepler orbit calculated for the center of mass. 



Accretion disks are composed of dust and gas and the dust tongas mass ratio is estimated 
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, l3| to range from ~ 10 2 to ~ 1. UV radiation photo-ionizes [J, [5], Q] the dust and gas so 
the accretion d,s k can be considered as a dusty plasma fl consisting of charged dnst grarns, 
electrons, and ions. The magnetorotational instability (MRI) [8] and the unipolar induction 
dynamo (UID) [9j] assume accretion disks are axisymmetric ideal MHD plasmas and neglect 
dusty plasma physics effects. The MRI and UID additionally assume that accretion disks 
obey both Kepler dynamics and ideal MHD. Thus accretion disks are considered to be 
ideal MHD plasmas in a circular Kepler orbit about a central object; they are supposed 
to conserve angular momentum and have frozen-in magnetic flux. One can then ask what 
does the trajectory of an individual particle in the accretion disk look like. Since MHD is 
assumed and MHD is based on the assumption that all particles make cyclotron orbits, it 
seems that the individual particle in question should be making a cyclotron orbit. On the 
other hand since the whole plasma is supposed to be making a circular Kepler orbit around 
the central object, then presumably this cyclotron-orbiting particle must also be making a 
Kepler orbit around the central object. A Kepler-orbiting cyclotron orbit is not an obvious 
concept to visualize, at least to this author. This conceptual difficulty suggests that instead 
of assuming that a particle is simultaneously Kepler-orbiting and cyclotron orbiting, one 
should go back to first principles to investigate how Kepler and cyclotron orbits relate to 
each other. Perhaps it will then become obvious how to visualize a Kepler-cyclotron orbiting 
particle, or perhaps not. 
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II. KEPLER EFFECTIVE POTENTIAL AND ORBITS 



Let us begin with a brief review of Kepler orbits. A particle a with conserved angular 
momentum L a = m a r 2 d(j)/ dt has radial motion in the Kepler "effective" potential [h]] 

L 2 m a MG 

XKeplerir) = (1) 

2m a r z r 

where M is the mass of a central object. Particles with energies at min XKepieri^) have 
circular trajectories with r = L 2 /m 2 MG, velocity vk = a/ MG/r and constant angular 
velocity = a/MG/t 3 whereas particles with energy exceeding mm\ Kpplpr (r) have elliptical 
trajectories and variable angular velocity d(p/dt = L a /m a r 2 [10(. 



III. GUIDING CENTER ORBITS IN COMBINED MAGNETIC AND GRAVITA- 
TIONAL FIELDS 



In contrast to Kepler dynamics 10j, basic plasma theory , Q shows that the 

guiding center of a charged particle in combined magnetic and gravitational fields, but no 
electric field, drifts at the velocity 



rrir. 



v 



9a q a B 2 



g x B (2) 



where g is the gravitational acceleration. For a mass M central object, g = 
MGV (r 2 + z 2 ) and so the guiding center drift in the z = plane is 

MG q Q i 

v 9<J = — <p = v K (j) (3) 

r io ca lo co 

where u Cf7 = q a B z /m a is the cyclotron frequency. It is seen that \ g is smaller than the 
Kepler velocity vk by Q/ \w ca \ , an enormous ratio for electrons and ions since \oo ce \ and \u C i\ 
are many orders of magnitude larger than Q for typical field strengths. If all the particles 
move at a much slower velocity than the Kepler velocity, then how could the center of mass 
move at the Kepler velocity? Furthermore, the \ ga add up to give the azimuthal current 
J g = Y2 n vQcr v gv = pg x B/5 2 where p = '^2,n (T m (T . The gravitational drift current results 
mainly from hea^ particle motion Q and gravity is balanced by the magnetic face 
(i.e., 3 g x B = — pg [13]) rather than by centrifugal force which is insignificant for this 
example. The ideal MHD Ohm's law E + U x B = actually fails here, because the Hall 
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term J x B/ne in the zero-pressure, generalized [14J] Ohm's law E + UxB = JxB/ne 
cannot be dropped since J/ne nearly equals the center of mass velocity U. 

What sort of trajectory does an actual charged particle follow in an astrophysical sit- 
uation? Is it the Kepler orbit assumed in Refs. IS, 9| or the much slower gravitational 

n 

drift derived in Ref. [13(? We show here that even though charged particles in a strong 



magnetic field can rotate at vk as assumed in Refs. |8|, , the motion is not governed by 
EqJT] so charged particles (and hence a plasma) do not in general obey Kepler dynamics. 
This analysis leads to the realization that dust grains having a critical charge to mass ratio 
spiral in towards the central object and so could provide a gravitationally powered dynamo 
suitable for driving astrophysical jets. 



IV. HAMILTON! AN DESCRIPTION OF ORBITS IN COMBINED ELECTRO- 
MAGNETIC AND GRAVITATIONAL FIELDS 



We consider the axisymmetric charged particle Hamiltonian 
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2 2 2m„r 2 V 2vr v ' ' V 



+ q a V(r,z,t) - 



m a MG 
\/r 2 + z 2 



Here ip = 2irrA ( i ) is the poloidal flux and is related to the magnetic field by 



B 



1 dip* dA z ~ 1 dip* 

-T — + —Z, 



(4) 



(5) 



lixr dz dr 2irr dr 
V is the electrostatic potential, and we note that ideal MHD ultimately comes from ap- 
proximations based on EqJD and not Eq(TJ Because of axisymmetry the particle's canonical 
angular momentum 



= m a r — + —i/)(r, z, t) = L a + —ip(r, z, t) 
at 2tt 2tt 

is invariant [16|, [17| . Equation H] is equivalent to the equation of motion 

dv 

m a — = q a (E + v x B) + m CT g 
at 



(6) 



(7) 



with E = -VV - dA/dt. 

We examine solutions to EqJT] in the z = plane for various charge to mass ratios, a 
uniform magnetic field B = B z z, and two representative V(r) profiles. In order to see 



the connection between the sense of particle injection and the magnetic field direction, the 
coordinate system definition we use here is such that positive <fi is determined by the direction 
of injection of the particle, i.e., the particle always has initial positive d<p/dt by assumption. 
This definition means that B z could be positive or negative since the direction of the z axis 
is determined by the sense of particle injection (i.e., we are defining the z axis so that the 
particle is always injected in the counterclockwise direction) and not by the direction of B. 
Because B z can be positive or negative, u cu will have the usual signage for B z > but will 
have opposite polarity from convention if B z < 0. To see whether or not particles make 
Kepler orbits we track particles starting with the same kinetic energy and L a as a neutral 
particle undergoing an elliptical Kepler orbit, a is defined to be the radius of the initial 
location and the x axis lies in the direction from the central object to this location. 

L a is conserved if q a = 0, but if q a ^ then rather than L a is the conserved 



quantity 



16 



171 ]. -B^ = will be assumed (to be justified later). Time is normalized to 
the Kepler frequency of a neutral particle undergoing circular motion at r = a, i.e., to 
fio = a/MG/o 3 and distances are normalized to a. The dimensionless variables are then 
f = r/a, t = Q()t, v = v/Qqci, L = L a /m a aHlo, H = H/mVt^a 2 , V(f) = q a V (r) / 'm^Vt^a 2 , 
and using = B z nr 2 , 

f 2 . (8) 



m a a 2 Q \dr 2Q 
The dimensionless z = plane Hamiltonian is thus 

Rotation of a plasma in a magnetic field polarizes the plasma radially and the resulting 
V(r) corresponds to the voltage to which the plasma capacitor is charged 17J]. There is 
thus no natural V(f) profile and hence no natural rotational velocity (e.g., see Refs. 171. l20|). 
We consider two representative cases, namely (i) V(r) = and (ii) V(r) = 2f 1 / 2 Lo> ccr /f2 . 
Case (i) corresponds to Eqf2] while case (ii) corresponds to the 'Kepler' equilibria assumed 



in Refs. 



HE 



9|. A third possibility, not discussed here (see Ref.[2l| instead), sets V(f) to give 
a rotation velocity equal to that of the central object (so-called "co-rotational velocity"). 

We consider all possible values of u ca /Q , namely |c<j ccr /fio| << 1) l^ccr/^o] >> 1, and 
|co> ccr /fio| = 0(1) with u ca /Qo either positive or negative. |co> C(7 /fio| >> 1 is typical for elec- 
trons, ions, and large charge to mass ratio dust grains, whereas |co> co -/f2 | << 1 corresponds 
to dust grains with very small charge to mass ratios or macroscopic charged particles such as 



spacecraft [22( . The voltage Vd to which a dust grain becomes charged depends on the charg- 
ing mechanism and the dust grain size; Vd typically lies in the range 1 volt < | Vd\ < 100 volts. 
Since the dust grain charge is Qd = 4:ne r d Vd, the charge to mass ratio Qd/rnd = ^oVd/r 2 d p d 
lies in the range 10 -4 — 10 2 C/kg for typical dust grain radii 0.1 /im < < 10 /im and 
typical dust grain intrinsic mass density pd = 10 3 kg m~ 3 . The dust grain uo c d/Q ratio is 
thus 9-15 orders of magnitude smaller than that of an electron and 6-12 orders of magnitude 
smaller than that of an ion. 

The last two terms in Eqj9] can be written as a normalized effective potential 

™ = \&-&i +v v-\- <io) 



If Q a — then ui ca = 0, V = and \ ~ * XKepier in which case Kepler dynamics [lOj is 
retrieved. 

However, when uj C(7 ^ and V(f) is arbitrary, the dynamics is non-Keplerian and Eq. 
[TU1 has minima when 



PA 2 ( u„\ 2 ldV 



f 2 J \2f2o / f df f 3 
By using EqJHl EqUH can be recast as 

2 



0. (11) 



uj C(J dd> 1 dV 1 . , 

dr J \lodr r or r A 

so a particle with H equal to the effective potential minimum has an angular velocity 




d(f) U! ca / / U! ca \ z ldV 1 



dr 2^ n V V 2 fin / f df f 3 



• --T- • -i- (13) 



Figure 1 plots x(f) and numerically calculated x-y plane trajectories for a range of uj ca /Qo 
values and for the two V(f) cases. In all trajectory calculations the particle initial position 
is x — l,y = and the initial velocity is v x = 0.4, v y = 1 (i.e., particles start at the same 
position with the same initial velocity and the same initial mechanical angular momentum 
L a ). The trajectories are calculated from r = to Att (i.e., two circular Kepler orbit periods) 
and the energy H is shown as a dashed line in the effective potential plots. An f = 1 reference 
circle (dashed) is shown in the trajectory plots. 

Figure 1 shows that the trajectory depends strongly on both the V(r) profile and on 
oj C(7 /VLq. The u ca /fl = situation (fifth row of Fig. 1) is a classic elliptical Kepler orbit 
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as prescribed by EqJT] and is independent of V(f) because a neutral particle is insensitive 
to electromagnetic fields. However, when uj cct /VLq is finite, Fig. 1 shows that the effective 
potential and trajectories differ qualitatively from the classic neutral particle effective poten- 
tial and elliptical Kepler orbit. It is therefore incorrect to characterize a plasma composed 
of particles with |u; co ./f2o| >> f as being in a Kepler orbit (as done in the MRI and UID 
models) because particles in such a plasma are not governed by Eq{T] and, for example, do 
not make elliptical orbits with the central object at one focus of the ellipse (such orbits are 
the "hallmark" of Kepler dynamics). 

Insight into the orbits shown in Fig. 1 can be obtained by examining solutions to EqJT3l 
For \u C(J /2VLq\ « I and V{f) = 0, Eq[J3]has the roots 

^ =± J (u) 

dr f 3 / 2 2Q 

so heavy charged dust grains with H equal to the minimum of the effective potential make 
circular orbits with a small uj cct /2VLq correction to the Kepler frequency. Heavy charged 
dust grains with H slightly above this minimum will make precessing elliptical Kepler orbits 
having small tu C(J /2f2o corrections (see uj C(J /VLq = ±0.1 cases in Fig. 1) and these corrections 
will increase with the charge to mass ratio. 

For |u; co ./2fio| >> 1 and V(f) = the two roots of EqJT31 are 



dcj) oo ca dcj) 1 VLq 
dr Qq ' dr f 3 u, 



(15) 



ca 



The first root corresponds to a so-called axis-encircling cyclotron orbit this root is not 
likely to be physically realizable in astrophysical situations since its corresponding azimuthal 
velocity exceeds the Kepler velocity by the large ratio |kj CCT /fi |- Normal cyclotron orbits 



correspond to a particle oscillating ly, ll7| in r about a local minimum of x{ r ) an d are 
associated with the second root in EqfT5l The second root is just v 9 prescribed by Eqj2] 
[ill . H . Q and, as discussed above, is smaller than vk by the factor |fi /^co-|- In reality 
l^o/^co-l would be so enormous that electrons and ions would have negligible azimuthal 
displacement during one Kepler period of a neutral particle. These slow drift orbits are 
shown in the top and bottom rows of case (i) in Fig.l. In accordance with Eqj2] heavy 
particles drift faster, negative and positive particles drift in opposite directions, and the 
drift velocity decreases as B z increases. 
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For case (ii), V{r) = 2r l / 2 u ca /VL, and so EqJT2l becomes 



dr 



^3/2 



dr r 3 / 2 



+ 



0J r 



^0 







(16) 



where one root is the circular Kepler-like orbit dip/dr = 1 with f = 1. Although the dtfi/dr = 
1 root looks superficially like a neutral particle Kepler orbit, the orbits are not elliptical, but 
nearly circular and, as in tokamaks and spheromaks, stay within a poloidal Larmor orbit of 



a constant ip surface 
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231 ] . The effective potential minimum has the same radial location 



as Eq{TJ but the profile is an extremely narrow trough with a large vertical offset (positive 
or negative, depending on the charge polarity), not a shallow broad well as for Eq{TJ 



V. ORBITS OF PARTICLES WITH ZERO CANONICAL ANGULAR MOMEN- 
TUM: DYNAMO FOR DRIVING ASTROPHYSICAL JETS 

A strange behavior is evident in the u cr7 /Qo = —2.0 row of Fig. 1: the effective potential 
goes to minus infinity on the left and the particle spirals inwards toward the origin in the 
x-y trajectory plots. This corresponds to P^ = and is unlikely for electrons or ions because 
they typically have |a^ co -/O 1 >> 1. However, P^ = could occur for dust grains because, 
being heavy, dust grains have \ui ca /2Q \ many orders of magnitude smaller than electrons 
or ions. As seen from EqJS], P$ = occurs if difi/dr = —uj C(7 /2Qo or, in un-normalized 
quantities P^ = occurs when dcj)/dt = —u Cf7 /2 in which case d<fi/dt also becomes invariant. 
In this situation the normalized effective potential, Eq JTUl reduces to 

* (f) 4(^) 2f2+ ^>-F (17) 

which has the remarkable feature that no centrifugal repulsion exists so the particle always 
falls towards f = with constant d<f)/dr (i.e., it spirals in). This differs qualitatively from 
P<t> 7^ particles which are constrained to orbit at a fixed average radius. The normalized 
radial force acting on a P^ = particle is 

dx 1 / uj ca \ 2 _ dV 



F = _r*=-_ ^ f- — -— (18) 

dr 4 V / dr f 2 K ' 

which is negative for any potential having dV/df > 0. Thus F is negative for both the 
V = and the 'Kepler' potential V(f) = 2f l ^ 2 uj ca /Vt. If there is a distribution of dust grain 
sizes, then a corresponding distribution of u C(T /Q values will develop and some subset will 



have uj ca /flo = — 2. The situation u ca /Qo = —2 could occur for negatively charged dust 
grains injected with B z > or for positively charged dust grains with B z < 0. This latter 
case is the likely one and if one were to adopt the common convention that the z axis is 
defined by the direction of B , this would correspond to retrograde injection. 

Dust grains accreting to a circumstellar disk will typically absorb stellar UV photons 4j 
and become positively charged by emitting photo-electrons [6|, |7|]. Some of these dust grains 
will satisfy d(f)/dt = —uj cct /2 if B z < and so have P<f> = P^ = 0. These dust grains have 
mechanical angular momentum L a = m a r 2 d(f)/dt at the instant before becoming charged 
by photo-emitting electrons, i.e., they have mechanical angular momentum L a = —q fT ijj/27r 
at the instant before they become charged. Since neither r nor d<p/dt is changed at the 
instant of charging, their canonical angular momentum P^ = L a + qip/27r becomes zero at 
the instant after charging. 

The infalling P^ = dust grains will accumulate at small f and create a positive space 
charge there. The photo-emitted electrons, stranded at large f (since electrons have P^ ^ 0), 
will create a corresponding negative space charge at large f. The positive and negative space 
charges will tend to cancel any polarization charge, e.g., the polarization charge associated 
with an initial V(f) = 2f 1 ^ 2 uj ca /Vl potential. Accumulation of infalling P<p = positive dust 
grains will eventually create an outward radial electric field E* (i.e., opposite direction to that 
associated with the V(f) = 2f 1 l 2 uj ca /VL potential). This accumulation will cease when E* 
becomes sufficiently large to create a force q a E* which cancels F. In un-normalized variables 
this cancellation occurs when E* = dV ot her/dr + {roj 2 a /A + MG/r 2 ) m a /q a where V Q ther is 
the potential profile that would exist due to particles other than the accumulating P^ = 
dust grains. The inward falling P^ = dust grains constitute a radially inward conduction 
current so J ■ E is negative and the system converts gravitational potential energy into 
available electrical power; i.e., it is a dynamo. The positive voltage near r = will drive 
bipolar axial electric currents / outwards from the z = plane. These currents will deplete 
the positive space charge which will result in a net force F — q^E* that will drive additional 
Ptf, = dust grains towards r = where they will replenish the positive space charge. Thus, 
the system continuously converts the gravitational potential energy of the Pa = accreting 



dust grains into a battery-like electrostatic potential which drives the po 



oidal current of an 



astrophysical jet. The jet itself is accelerated by the dB?/dz force [2J, |25|, |26j associated with 
the axial non-uniformity of the jet poloidal current I(r,z) since B,p(r,z) = fioI(r, z)/2itr. 
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The anti-symmetry of / with respect to z means that = at z = consistent with the 
assumption made earlier. 

We note that the ability of = particles to cross magnetic flux surfaces has recently 
been observed in a laboratory experiment (27) ] . 

The assumption used in this paper that B z is spatially uniform is a simplifying idealization 
that enables the analysis to be both brief and focussed on the distinction between Kepler 
and charged particle orbits. However, an actual accretion disk will almost certainly have 
B z depend on both r and on z so in order to satisfy V • B = there will also have to be a 
B r (r, z). This indicates that the axisymmetric magnetic field would be best described using 
a poloidal flux function ip(r, z) i.e., B(r, z) = (27r) _1 (Vt/> x V</> + HoI(r, z)V(fi) . This more 
general description of the magnetic field has been used in Ref.|28], a much lengthier analysis, 
where three dimensional particle orbits in an approximately self-consistent magnetic field are 
considered using a generalization of the Hamiltonian method presented here. Specifically, 
the poloidal flux function ip(r, z) in Ref. 28( results from a toroidal current due to toroidal 
motion of charged particles and particles are not restricted to the z = plane as in the 
present paper. 

The dust grains might be so densely packed as to be optically thick in which case photons 
from the central object would not reach the dust and the dust would not become charged; 
this would constitute a so-called 'dead-zone'. The condition for a dust cloud to be optically 
thick is naL > 1 where n is the dust density, a is the dust cross-section, and L is the 
characteristic length scale of the dust cloud. However, the condition for a dust cloud to 
be collisional is also naL > 1 and so dust grains in an optically thick cloud would be 
collisional. This collisional, optically thick state would likely be transient, because collisions 
are expected to cause coagulation of the dust grains [lj in which case their radius will 
increase. Since the mass of an individual dust grain is rrid = 4'np d r 3 l /3 and since coagulation 
does not change the total mass M of all the dust grains, the number N = M/m of dust 
grains and hence the density n ~ N/L 3 of dust grains scales as r^" 3 . Because the dust grain 
cross-section a scales as r^, the product na scales as r^ 1 and so na decreases as a result of 
coagulation. Coagulation of dust grains will thus reduce na until naL becomes less than 
unity in which case the dust cloud will become collisionless and optically thin. At this point 
dust would become charged (ionized) via photo emission and commence the collisionless 
trajectories discussed here. The dead zones would thus disappear as a result of coagulation. 
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This issue is discussed in more detail in Ref. 

Ref.|28| discusses several other issues including: the charging rate of dust grains (i.e., 



effective ionization rate of dust grains), collisions of dust grains with gas and other dust 
grains, and the topological properties of astrophysical jets. These various issues are used to 
define a parameter space for astrophysical jets powered by the gravitational energy released 
by accreting = dust grains. A self-consistent set of parameters is given for the example 
of the circumstellar accretion disk of a young stellar object (YSO). 
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FIG. 1: Effective potentials and x-y plane trajectories for particles starting at x = 1, y = with 
initial velocity v x = 0.4, v y = 1 with range of u) C(J /VLq values and two V(f) cases. Elliptical Kepler 
orbits (i.e., Eq{T] effective potential) occur only when Lo ca jVt = 0. Particles with u) CC7 /VIq = —2 fall 
towards f = and have = 0. No particles malle Kepler-like elliptical orbits when |cj ccr /r2o| >> 1- 



